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1. Introduction 

In the past twenty years we have witnessed steady increase in the interest in gravitational 
waves from astrophysical sources. Specially driven by the design and construction 
of laser interferometric detectors, both ground and space based. Alongside with this 
experimental developments theoretical progress has also been steady. We are now 
in conditions to predict the gravitational radiation from the astrophysical scenarios 
expected to produce the strongest signals, i.e. the merging binary black holes. 

The two main scenarios involving black holes are, first, galactic binaries with black 
holes having comparable masses (a few solar masses). They are, for instance, the product 
of a supernova explosion plus a subsequent accretion. The second scenario involves a 
supermassive black holes (with several million solar masses) residing in the center of 
an active galaxy. They attract stars in the inner nuclei towards unstable orbits with a 
subsequent plunge generating observable gravitational radiation. This scenario clearly 
involves extreme mass ratio collisions. Let us also mention that a less common event, 
but most energetic, is the close collision of galaxies and consequently of supermassive 
black holes in their respective cores. 

From the theoretical point of view one advantage of dealing with binary black 
holes is that one can treat the problem of generation of radiation in terms of only its 
gravitational field, ignoring the (small) effects of matter around the binary system. A 
second important feature is that the equations of General Relativity scale with the 
total mass of the system. In this way, one can choose the dimensions of the systems a 
posteriori, i.e. after solving for the scale free problem (See |T] for the case of three black 
holes.). It is then convenient to characterize the binary black hole systems in terms of 
the mass ratio of its components (besides the individual spins, orbital parameters and 
spatial orientation with respect to the observer). 

In binary black hole systems most of the generation of gravitational radiation take 
place during the final few orbits before the merger. In the case of comparable masses, this 
stage involves highly nonlinear interactions among black holes and can only be described 
by directly solving numerically the full General Relativity field equations. Until recently 
this represented an insurmountable task that hold the field for nearly thirty years, but 
during the year 2005 two successful approaches [21 [3j Hj have lead to stable codes that 
allow to simulate binary black holes in supercomputers. This breakthrough in Numerical 
Relativity lead to numerous studies during the last year, including the last few quasi- 
circular orbits [3 [6], the effect of ellipticity [7J, and spin-orbit coupling leading to a 
change in the merger time [Sj, corotation [9] and spin- flip and precession [10]. Finally, 
unequal mass black holes have been studied in [TU [T21 [131 Ej reaching a minimum mass 
ratio of nearly 1:4. But simulations with mass ratios up to 1:10 are currently underway. 
Generic binary simulations, i.e. unequal masses and spins, have first been reported in 
[15], this simulation lead to the shocking discovery that merging spinning black holes 
can acquire recoil velocities up tp 4000 km/s [16]. Even multi-black hole spacetimes are 
now possible to evolve numerically [T7J, [18] . 
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In the small mass ratio regime the smaller hole orbiting the larger is considered 
as a perturbation. In this approach the smaller hole is described by a Dirac's delta 
particle and the spacetime is no longer empty but has a non- vanishing energy-momentum 
tensor at the location of the particle. The simplicity of this approach is appealing, 
but the problem notably complicates when self-force effects are taken into account to 
compute the correction to the background geodesic motion. A consistent approach to 
this problem has been first laid down ten years ago by Mino, Sasaki and Tanaka [19] 
and later confirmed by Quinn and Wald [20j by providing a regularization procedure. 
(See [2TI [22] for a detailed review and [23j for a practical regularization method.) 
Self-force corrections can be considered second order effects on the mass ratio of the 
holes, the natural perturbation parameter for the system. To consistently compute 
the gravitational waveform, radiation energy, angular and linear momentum radiated 
to infinity (and onto the horizon) we need to proceed to solve the second order 
perturbations problem for the gravitational field. The formalism to study second (and 
higher) order perturbations of rotating black holes with sources was extensively discussed 
in [23] based on the Newman- Penrose approach of curvature perturbations. 

The first explicit computation of the gravitational self-force was done for the headon 
collision of two black holes in Ref . [25l [26] . This allow us to complete here, for the first 
time, the second order program applied to the headon collision of black holes. Note that 
this step is also crucial in order to close the gap between the full numerical simulation 
with comparable masses and the perturbative approach for extreme mass ratios. 

Here, the second order calculation is required to derive precise gravitational 
waveforms which are used as templates for gravitational wave data analysis. In general, 
this second order computation has to be done by numerical integration. It is hence 
important to derive a well-behaved second order effective source, and we will focus on 
this problem in this paper. The second order analysis was pioneered by Tomita [27| 28J, 
and vacuum perturbations in the Schwarzschild background was studied by Gleiser et 
at [291 SOI [31], [32] . There are also studies on second order quasi- normal modes [33] [33], 
and the second order analysis was also extended to cosmology [35] [36] [37] [38] [39] . 

This paper is organized as follows. We consider second order metric perturbations 
and the equations they satisfy, i.e., the perturbed Hilbert-Einstein equations of General 
Relativity in section [2] In section [3] we discuss the first order metric perturbations in 
the case of a particle falling radially into a Schwarzschild black hole by using the Regge- 
Wheeler-Zerilli formalism [401 |4"T] . This formalism in the time domain is summarized in 



Appendix A[ To calculate the second order source, the first order I = mode is given 



in an appropriate gauge in Appendix C In section [4] we derive the regularized second 
order effective source where the singular behaviors are removed analytically. In Sec. [5] 
we summarize this paper and discuss some remaining open problems. Details of the 
calculations are given in the appendices. Throughout this paper, we use units in which 
c = G = l. 
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2. Second order metric perturbations 

We consider second order metric perturbations on black hole backgrounds, 

g»v = g^v + h$ + hf v , (l) 

with expansion parameter fi/M corresponding to the mass ratio of the holes. Here, g^ u 
is the background metric, and superscripts (i) (i = 1, 2) denote the perturbative order, 
i.e., h$ and are called the first and second order metric perturbations, respectively. 
In the perturbative calculation, we raise and lower all tensor indices with the background 
metric. The Hilbert-Einstein tensor up to the second perturbative order is given by 

GAU = + + hF>\ , (2) 

where we have omitted the spacetime indices /i and v of the metric perturbations, 
h$ and Jiffi, and ignored 0((/z/M) 3 ) and higher order terms. G$ is the well known 
linearized Hilbert-Einstein tensor, 

G^l[H] = —-H tlu . a ' a + H a ( KU y a — R a ^/3 U H al3 — -H a a ]lJiU — -gf JtlU (H\a aX — H a a ; \' x ) , (3) 

Here, H^ v denotes Kju, or h$, and semicolon ";" in the index denotes the covariant 
derivative with respect to the background metric. G^l consists of quadratic terms in 
the first order perturbations, 

G«[/i« *«] = B$[hP\hF>\ - i v 4 2)Q [^ (1 U (1) ] ; (4) 
n$[hVhV) = \h%,»h^ w + \h^{hX, + C - 2 <U) 

On the other hand, the stress-energy tensor includes three parts; 

T _ 1 T (2,SF) , T (2,h) (c\ 

The first order stress-energy tensor T$ which is the one of a point particle moving 
along a background geodesic, is given by 

/ + OO J fl J V 

S^^ x -z{r))^-^-dr, (6) 
dr dr 

where 

z» = {T(t),R(t),Q(t)Mt)}, (7) 
for the particle's orbit. T^V denotes the deviation from the geodesic by the self-force 



as derived by the MiSaTaQuWa formalism [19l [20] . We do not treat this stress-energy 
tensor explicitly in this paper. And T^' h \ which is purely affected by the first order 
metric perturbations, is written as 

= -\n rh^6^(x-z(r)) d -^ d £dr, (8) 
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where we have used the determinant 

g = g{l + h^ a ) , (9) 

up to the first perturbative order. 

With the above expansion of the Hilbert-Einstein and stress-energy tensors, we may 
solve the following equation for the first perturbative order, 

= 8 7tT«. (10) 

And for the second perturbative order, we have the following equation. 

^1 = 8, (^ +3 y))-Gg[fc( 1 U«]. (11) 

Once the first order metric perturbations h^- 1 ' (and the self-force) are obtained, we may 
solve (Till with a second order source that can be considered as an effective stress- 
energy tensor. Systematically expanding the Hilbert-Einstein equations, one can obtain 
the perturbative equations order by order |32J EH H3J 0U H5] . 

3. First order perturbations in the Regge-Wheeler gauge 

In this paper, we consider the Schwarzschild background, 

ds 2 = -(l -^pjdt 2 + (l 1 dr 2 + r 2 (d6 2 + sin 2 dd<p 2 ) . (12) 

in Boyer-Lindquist coordinates. 

Before considering the second perturbative order, it is necessary to discuss the first 
order metric perturbations, i.e., the first order Hilbert-Einstein equation given in (Tl0|) . 
The Regge-Wheeler- Zerilli formalism [4"0],|4"T] is used here. The basic formalism has been 
given in Zerilli's paper [41], and it has been summarized in the time domain in [4H1 Wf\ . 



In |Appendix A[ we establish our notation and summarize the Regge-Wheeler-Zerilli 



(i) 

fj,v 



formalism in the time domain. 

The treatment of the first perturbative order is as follows. First, we expand h 
and T$ in ten tensor harmonics components, given in (lA.ljl and (1A.2I) . We then obtain 
the linearized field equations for each harmonic mode. Here, for example, for the even 
parity modes which have the even parity behavior, (— l) e under the transformation 
(9, (f>) —>■ (n — 9, cf) + 7r), we may consider the Zerilli equation in (|A.8[) . Finally, imposing 
the Regge-Wheeler (RW) gauge conditions: 

h (e)(l)RW _ h (e)(l)RW _ r (l)RW _ n „\ 

a oem — 'hem —^im — u > \ LO ) 
where the suffix RW stands for the RW gauge, we obtain the first order metric 
perturbations as in (IA.110 . 

3.1. Geodesic motion and the first order stress- energy tensor 

We consider a particle falling radially into a Schwarzschild black hole as the first order 
source. Assuming O(r) = $(r) = 0, the equation of motion of the particle is given as 



f dR(t) \ 2 ( 2M\ 3 1 / 2M\ 

\nr) V W)J E 2 + V W)J 



TS+ l--57^ . ( 14 ) 
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where R(t) is the location of the particle and the energy E is written by 

/ 2M\ dT(r) , s 

We will also use 

d 2 R(t) __ 3 / 2M \ 2 M 2 / 2M \ M 



dt 2 E 2 y R(t) J R(t) 2 + 2 V 1 R(t) ) R(t) 2 ' (16) 

to simplify equations. The tensor harmonics coefficients of the first order stress-energy 
tensor which are given in table IA1| 



= „ 1 ^§ M ^^S(r - R(t))Y L (0, 0) , 

41(t, r) = V2i» R ^ m 2M dF f ±6(r - R(t))Y L (0, 0) . (17) 

The remaining coefficients are zero. Because of the symmetry of the problem, we have 
only to consider the even parity modes. 

3.2. First order metric perturbations (£>2) 

We introduce the following wave-function for the even parity {£ > 2) modes, 
i>Z n (t,r) = ^^[KjST '(t,r) 

(r-2M) f T Ai)RW ( . \ d (i)rw. 



This function ipZ n obeys the Zerilli equation, 

■^ + ^T2- Vr n (r)) r) = SZ n (t, r) , (19) 



where r* = r + 2Mlog(r/2M — 1), the potential Vf ven is given in (1A.9j) . The source 



S|™ n is calculated as in ( 1A.10I) . For the £ = and 1 modes, we will need a different 
treatment as described in section [3^1 (See also [151 B9] .) 

The reconstruction of the first order metric perturbations from this wave-function 
in the RW gauge have been given in (lA.llj) . In the head on collision case, the metric 
perturbations in the RW gauge are C° (continuous across the particle). One can see 
this as follows. First, using the following linearized Hilbert-Einstein equations for each 
harmonics mode, 

rr(l) _ rr(l) o„2r(l) 

n Qlm n 2lm _ Q/ " s hn (nQ) 

2 " v /£(£+l)(£-l)(£ + 2)/2' 

and = 0, we obtain H^)^ w = H^)^ w . Then, removing Hf^ W from two linearized 
Hilbert-Einstein equations , 
d_ 

Or 



1 — ) H \(n, 



(l)RW 



9 ,„(l)HW . jf(l)RWs »(1) 



_(irVJ nrv i zr\ L ) nvv \ — Ky L > (o-\\ 
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(1)RW 

im 



dt 
1 



+ 1 



2M 



dr 



i u { 



im 



K 



Im ) 



2M 



+- 



(1)RW 



-H, 



8ti (r 



+ 

r 

2M) 



2 "Qfm 



we obtain the relation 



2M 



3^ 
<9r 2 



(l)iW 
2 fen 



V^ + l)/2 



(i) 



(22) 



(t,r) 



<% 2 



+ 1 



2M 



<9r 2 



+ (lst differential terms of ffj^ an d ^im^' 



(23) 



where we have used £>, 



(i) 

O&ra 



B 



(i) 



0. Therefore, we find that H^^ W and Kl} RW have 



the same differential behavior. Here, we note that the wave-function ipg^ n behaves as a 
step function around the particle location because of (1191) (and (j24p ). Thus, it is found 



(1)RW 

£m 



6{r — R(t)) with use of (1181) . This means that K^ RW is C°. From 



that d r K L 

O with the above fact, d r H { $* W ~ 6(r - R(t)) is derived , i.e. H^ W is also C°. 
(See [25].) We note that we can take up to second derivatives of the function with 
respect to t and r around the particle location as in (1271) . 

In the next subsection, we treat only the I = 2 mode which is the leading 
contribution in the first order perturbations. In this 1 = 1 mode, we may consider 
only the m = mode because of Yi m (0, 0) = for m ^ in ( ITTj) . 



3.3. i = 2, m = mode 



We focus here only on the 
(fl9j) becomes 



2, m = mode. For this mode, the Zerilli equation in 



<9 2 



<9 2 



<9t 2 + dr* 2 



-6 



- 2 M) (4 r 3 + 4 r 2 M + 6 rM 2 + 3 M 3 ) 
r 4 (2r + 3M) 2 



i 7T 



// (2 i?(t) 2 - 2 R(t)E 2 M + 6 i?(t)M + M 2 ) (i?(t) - 2 M)' 



Er 3 (2r + 3M) Z 



xK,* (0,0)5 (r-R(t)) + — 



8 7T ii{R{t)-2MY 



3 ER{ty {2R{t) + 3M) 



Y 2 * (0,0) ^-5 (r-R(t)) , (24) 



where we have used the formula 

F(r)4-5'(r - i?) = F(R)^-5(r - R) - ^-F'ir) 
dr dr dr 

to simplify the source term. 

Here, we decompose the wave-function in the following form, 

9%*{t, r)6(r - R(t)) + *%(t, r)6(R(t) - 



dr 



5{r - R) 



(25) 



9»(t,r)0{r-R{t))+9»(t 



H i 



*»(t,r) = r)-*£(t,r), ^ (t, r) = r 



(26) 



where vE^o* anc ^ ^20 ' i- e -' s ^ so ^ & anc ^ are homogeneous solution to the Zerilli 
equation. Using the fact that the first order metric perturbations are C , the following 
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six quantities can be derived 
« / 8tt uE R(t) „ . 

<9r 20V ;|r KW (2i?(t)+3M) 2 (2M-fl(t)) 20 V ; 

d e fiE R{t) dR(t) , 

-* 20 (t,r)| r=i?w = -Stt (2 - 3 M) (2 - _ ^ — y M (0, 0) , 

9 2 fiE (7 M 3 — 4 i?(t)M 2 + 4 R{tfM - 8 i?(t) 3 ) 
— % (t,r)\ r=m - -Svr (2it > (t) + 3 M) 3 (2M-^)) 2 ^ ( °' ° } 

:^20l^ r )|r= J R(t) = 87T „ ——2 — 1 20 V^' , 



^ *n w (2i?(t) + 3M) 2 (2M- J R(t)) 2 d* 

4?*I>fo | r= * w = -8tt y* (Q Q) (27) 

dt 2 20K Jlr H(t) (2R(t) + 3M) 3 R(ty 20 K ' K ' 

In the above equations, first we take the derivatives, and then set r = R(t). These 
quantities allow us to calculate the coefficients of the 5 terms in the second order source. 

3.4- £ = mode (Monopole perturbation) 

Next, we consider the £ = perturbation (£ = 1 modes can be completely eliminated in 
the center of mass coordinate system) which is present only in even parity. The metric 
perturbations and the gauge transformation are given as 

h$o = (l - #ooo(*: r ) a ooo - iV2H$ (t, r)a 100 

1 - — J H^ (t, r)a 00 + V2K&\t, r)g 00 , (28) 

Co = {V (l \t,r)Y 00 (9, ( f>),V} 1 \t,r)Y 00 (9,<f ) ),0,0}, (29) 

respectively The metric perturbations transforms under the above gauge transformation 
from the G gauge to the G' gauge as 

H™\t, r) = H^(t, r) + 2^V™^(t,r) + -^-^(^(t, r) , (30) 
H%°\t, r) = HH%, r) - _L_ » (*, r) + (*. , (31) 

Sffifftr) = fl^(t,r) - + -^L_V^ (r) , (32) 

A-<; )G '(«, r) = /C G (i,r) - V« 1,G - G '(«, r) . (33) 

Here, we can choose Vq 1 ^ G ~^ G and V^ G ~* G so that H^qq = K^ z = where 
the suffix Z stands for the Zerilli gauge [JT]. In this gauge, the two independent field 
equations are given by 

dH ™&r) 1 (1)z _ 8vrr 3 (1) 

0r + r - 2ikT 200 ^ (r - 2M)^ ooott ' rj ' ld4j 
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d4l Z (t,r) 



1 



H^ Z (t,r) = -8n r A^(t, r 



(!)/ 



dr r- 2M 



(35) 



where A q anc ^ ^-oo are gi ven i n (HTI) . We solve the first equation, and then obtain 



r-2M 



Y* (0,0)6(r-R(t)). 



(36) 



Next, substituting the above quantity into the second equation, we obtain 



H^ Z (t,r) = 8^E 



1 



_ 1 R(t) 2 ( dR(t) \ 2 

r-2M R(t) - 2M (R(t) - 2Mf \ dt J 



xY* (0,0)6(r-R(t)). (37) 

It is difficult however to construct the second order source from the above metric 
perturbations, since these are not C°. We instead consider a new (singular) gauge 
transformation, chosen to make the metric perturbations C°. 

We consider the following gauge transformation. We call this the C gauge, where 
the first order metric perturbations is C° at the particle location 

2 7T fi Y* (0, 0) (r — 2 M) (r + 2 M) (r 2 + AM 2 ). 



3E 



2nfiEY* (0, 0) 



(r-R(t))R(tf 



dR(t) 



-INT(t) 
-r 2 R(tf - R(t) 2 r 3 - r A R(t) 



3 (R(t) - 2 M)\\r - 2 M) dt 

-10 r 2 MR{tf + 10 MR(t)r 3 + 8 Mr 4 + 10 rMR{tf - 16 rM 2 R{tf - 16 M 2 R{t)r 2 

(38) 

(39) 



•16 M 2 r 3 - R(tfr + 10 R{tf - 62 R{tfM + 72 R(t) 3 M 2 ) 6(r - R(t)) , 



{t,r)=4n l iEYJ a (a, 0) 



(r-2M) (r-R{t))R(tf 



INT(t) 



r 6 -2M) 2 
(-54 M 2 + 39 i?(t)M - 50 ME 2 R{t) + 12 E 2 R(t) 2 



9{r - i?(t)) ; 



-6 R{t) 2 + 16 M 2 E 2 ) I (R(t) - 2 M) 3 

M 

R(t) 



dt 



(-1 + 2E 2 )E 
6 , — arctan 



^fT^E 2 



+12 E In 



AM 2 
1-E 2 

AME 



+ 



1-E 2 
2R(t)M 



1-E 2 



— AR(t)M 



y/T^E 2 



R(t) 



2M 
1 - E 2 



-R(t) 



R(t) 



■,1/2 



2M 
1-E 2 



-R(t) 



1-1/2' 



(M(R(t) -2M)Y l 



(13 -R(t) 2 + 48 M 2 — 56 R(t)M) E[E 2 R(t) — R(t) + 2 M] 1//2 y^R(t) 
+ (i?(t) -2M) 3 ' ^ ^ 

We then obtain the £ = mode of the metric perturbations in this C gauge by using 
(EQD, (EH), (J325 and (J33J), 

4 7r /i F * (0, 0) (r — 2 M) (r + 2 M) (r 2 + 4 M 2 ^ 



E (R{t) - 2 M) 3 r 4 



-5AM 2 + 39 R(t)M 



-50 ME R(t) + 12 E 2 R(t) 2 - 6 i?(t) 2 + 16 ME) 6(R(t) - r 
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-792 R(tfr 3 M 3 + 792 R(t) 2 M 3 r 4 



4 7T/xF * (0, 0) 1 700 DW3J W 3 , 7no D/ri2» # 3J 



3 (i?(f)-2M) 2 r 7 (r-2M)£ 
+192 R(t)r 3 M 5 - 864 r 3 M 6 + 6 R(t) 7 MrE 2 - 398 R(tfr 3 M + 256 r 3 M 6 E 2 
-128 r 4 M 5 E 2 - 6 R(tfMr 2 E 2 + 308 R(tfMr 3 E 2 - 468 R(tfr 3 M 2 E 2 
-72 R(t) 3 r 3 M 3 E 2 - 144 R(t) 2 r 3 M A E 2 - 672 R(t)r 3 M 5 E 2 + 336 R(t)r 4 M 4 E 2 
+50 i?(t) 6 r 3 - 44 #(t) 5 r 4 - 24 r 7 M 2 E 2 - 50 R(tfE 2 r 3 + 72 R(tfr i M 3 E 2 
+24 R(t)r 7 ME 2 + 476 R(t) 3 r 4 M 2 E 2 + 12 R(tfrM 2 E 2 - 282 R(tfr 4 ME 2 
-6 R(t) 2 r 7 E 2 - 96 R(t)M 4 r 4 + 359 R(t) 4 Mr A - 968 R(t) 3 M 2 r A + 432 r 4 M 5 
+ 1022 R(tfr 3 M 2 - 12 R(t) 7 M 2 E 2 + 44 R(tfr A E 2 ) 9(r - R(t)) , 
^ = 128,^o(0.0) M'(r-2M) 



4.^(0,0) (r -*(«))*(*)■ rf fl(()( _ 72MV -288fl( ( )M= 



3 r 6 (r - 2 M) (i?(t) - 2 M) 4 dt' 

+68 rR{t)M 2 + 60 MR{t)r 2 - 132 MR(tfr + 144 rM 3 - 40 ,R(t) 3 M + 25 i?(t) 3 r 
+248 R{t) 2 M 2 - 12 i?(i) 2 r 2 ) 0(r - , 

ff<jL°(t,r) =87r/iSK^(0, 0) r - ^ 5 (-21 RitfrM + 22 #(t) 6 M 2 

; (r-2M)r 7 (i?(t) -2M) 2V 



+5 i?(t) 6 r 2 - 4 R{tfMr 2 + i?(t) 5 r 3 + 4 R{tfrM 2 + 4 R{t) 4 M 2 r 2 + i?(t) 4 r 4 

-4 R(t) 4 Mr 3 - 4 R(t) 3 r A M + 4 R(tfr 3 M 2 + ii!(i) 3 r 5 - 4 R(t) 2 Mr b 

+4 R(t) 2 r A M 2 + i?(t) 2 r 6 + 4 R(t)M 2 r 5 - 4 R(t)Mr 6 + 4 M 2 r 6 ) 0(r - , 

r) = -8 „ „ fiiatO. 0) - 2 fT y y » ^ - *M) . (41) 

Note that all of the above metric perturbations are C° at the particle location and go to 
zero at r = oo and the horizon. Using the metric perturbations in the C gauge, we have 
discussed the second order source in [50]. But, with this method, we can not obtain the 
second order gravitational wave at infinity because this is not an asymptotic flat gauge. 

We therefore must consider another treatment. The metric perturbations in (1411) 
have the following asymptotic behavior for large r 



H, 



oTfc r) = f 8 ^ + 16 ^) y *(0, 0) + O(r-) , flgjfa r) = 0{ r - A ) 



i^,r) =(^8 ^ + 16^^ (42) 
Hence, the metric up to the first perturbative order in the system becomes 
ds 2 = " (l " (1 " #ooo°(t, r)Y m (9, <P))dt 2 + r)Y 00 (6, <f>)dtdr 

+ 1 (1 + H ^(t, r)Y 00 (6, <j>))dr 2 + r 2 (l + K { £ c {t, r)Y 00 (6, 0))rffi 2 



r 

-l 



2M + 2^ rft2+ / 1 _2M + 2^ dr2 + rW . (43) 
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We thus find that this perturbation is related to the mass increase of the system. 

From the above analysis, we define the total mass M tot = M + /x E as that of the 
system. This means that the particle's mass is absorbed in the background Schwarzschild 
mass. Then, the first order £ = renormalized metric perturbations become 

4T(t, r) = 4l C (t, r) - H^ M (t, r) , 
HZ(t, r) = HH C (t, r) - H^ Q M (t, r) , 

H[l N (t, r) = H[l C (t, r) , K<£ N (t, r) = K$ C (t, r) , (44) 
where we have labeled this renormalized metric perturbations by N, and 

H$o M (t, r) = Hgjfit, r) = 8 ^r^Y* (0, 0) , (45) 

Now for the asymptotic behavior for large r, we have 

fl^(t,r) = 0(r-*), H^ N (t,r) = O(r^). (46) 

Next, we will use the coefficients of the first order metric perturbations labeled by N in 
to derive the second order source. 



4. Second order perturbations in the Regge- Wheeler gauge 



4-1. Second order Zerilli equation 

Since the first order metric perturbations contain only m = even parity modes, we 
can discuss the second order metric perturbations via the Zerilli equation. And we 
will choose the RW gauge condition. Here, we use a wave-function for the second 
perturbative order, 



1 



d 



r~-KT W (t,r) 



(r-2M)HH RW (t, 



(47) 



2r + 3M 

This is the same definition as in flA.7j) for the first perturbative order. Here, we have 
considered the contribution from the I = and 2 modes of the first perturbative order to 
the t = 2 mode of the second perturbative order since this gives the leading contribution 
to gravitational radiation. This Zerilli function satisfies the equation, 

r - 2 M)(4r 3 + 4r 2 M + 6rM 2 + 3M 3 )' 



o-even 



d 2 d 2 

I 6 

dt 2 dr* 2 



r 4 (2r + 3M) 2 



5 2 Z (t,r) 



— $20 (*> r 
8tt V3(r 



(48) 



2M) Z 9 o(2) 8n(r-2MY d (2) 



3(2 r + 3M) 
8V37r(r-2M) d 

3 di' 20 
8 \[2m (r — 2 M) (5 r — 3 M) M 

r (2r + 3M) 2 
32 y/Zm (3M 2 + r 2 ) 



2r - 31/ d t A2o{t}r) 



4V2 in (r — 2Mf d 



2r + 3M dr A ^ r) 
(2) r+ ^ _ 8 y / 3z7r(r-2M) 2 g (2) 
3(2r + 3M) <9r 020 



3r(2r + 3M)' 



2M) g(2) 



•^1 2o(^> r ) 
11 2o(*J r ) ■ 



(t,r) 
(49) 
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(2) 

The functions B 2 q etc. are derived from the effective stress-energy tensor on the left 



hand side of (fTTI) . 

T^f n) = M 2 /^ + TM>) - [*« ,/!«], (50) 

by the same tensor harmonics expansion as for the first perturbative order. 

The second order metric perturbations from X20 * n ^ ne g au S e are given by 

^A 20 (i,rj-b (2r + 3M)r2 Xaol*,rj+ ^ ^ol^j 

4^r(r-2M) (2) 8^3 ZTrr (r - 2 M) (2) 

+ 2r + 3M ■ A iao^' r J + 3(2r + 3M) ^2ol^j, l^J 

1^ (t, r) = r^*ff~ (*, r) + 3 (*, r) - ^^|x!o (*, r) 

8^ . ( 2 ) 



m^ (t,r)r, (52) 

(«, r) = Fiir r) + ^ w jg> (t , r) r2 . (53 ) 

Since the RW gauge is not asymptotically flat, we need to derive the second order metric 
perturbations in an asymptotic flat (AF) gauge to obtain the second order gravitational 



wave at spatial infinity. This is discussed in Appendix C 



Note that in (jSJ). the delta function 8^(x — z(r)) in T$ h ' includes an angular 
dependence, 5^(Vl — fi(r)) = J2e m ^m(^)^m(^( r )) • We h ave considered only the 
contribution from the £ = and 2 modes of the first order perturbations. Consistently, 
we must use only the three components with the factors, ha (£ = 2)F 2 m(^)^2m(^( r ))) 
f$ a {t = 2)Y 0m (n)Y^ m (Q(r)) and hg )a (£ = 0)Y 2m (Q)Y* m (n(r)) . 

In the second order source, given in ()49p . we may wander if there is any 5 2 term 
which prevents us from making the second order calculation. The answer is "No" . This 
is because, in the head-on collision case, the first order metric perturbations in the 
RW gauge are C°: G^[h^\h^] includes second derivatives and we need one more 
derivative to construct the second order source of (149j) . Here, (h^) 2 is C° x C° and 
its third derivative yields C° x 5' and 9 x 5 as the most singular terms. On the other 
hand, T^, includes only C° x 5 terms. Note also that there is no 5 2 terms coming from 
Tj& SF \ (which we ignored otherwise in this paper.) Using the result of [251 126]. we can 
include the contributions of T^l' SF \ 

In the following subsection, we derive the second order source of the Zerilli equation 
in ( 14*9]) . The summary is given here. From (j4"9"l) . we obtain the second order source as 

Sl^r) = ^SUt,r) + ^SUt,r)^ (54) 

where (2 ' 2) <S| and (0 ' 2) <S| are the contribution from (£ = 2) • (£ = 2) and (£ = 0)-(£ = 2), 
respectively. Note that while the above source term is locally integrable near the 
particle's location, some terms diverge as r — > 00 or 2M. This is not readily suitable 
for numerical calculations. We then consider some regularization for the asymptotic 
behavior. (See e.g. [31].) In order to obtain a second order source which behaves well 
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everywhere, we define a regularized Zerilli function by 

X 2 Z o(t, r) = xl(t, r) ~ X r 2 o g ' (2,2) " • (55) 

The best suited equation to solve the Zerilli equation numerically is then 

Zr n xl(t,r)=S^ ce (t,r), (56) 
where the regular source S^q 05 is given by 

S 2 Z reg (t,r) = (^S 2 Z (t,r) -ir n X 2 o S ' (2 ' 2) (t,r)) 

+ ((°' 2 >S 2 z (t, r) - Zr n xlT i0 ' 2 \t, r)) • (57) 

Regularized second order source from (£ = 2) • (1 = 2) 

When we consider the asymptotic behavior of the second order source for large r, we 
use the retarded solution of the first order wave-function ip^ n with the retarded time 
t — r*, which we expand in inverse powers of r. The wave-function becomes 

r) = F'j{t - r*) + -Fj(t - r*) + 0(r- 2 ) , (58) 
r 

where Fj is some function of (t — r*) and Fj(x) = dF(x)/dx. 

On the other hand, the wave-function is expanded near the horizon as, 

r) = F' H (t + r*) + ^F H (t + r*) + ^['^ F H (t + r*) 

+ 0((r-2M) 2 ). (59) 

Using these expansions , we derive a second order source which is regular everywhere. 

In order to obtain a well behaved source for large values of r, we define a 
regularization function by 

reg,(2,2) f s _ ^ / 9^l) M( A .(IJW, x ^ 

Note that the regularization function is not unique. Therefore, this affects the formal 
expression of the second order gravitational wave as in flC27[) . However, for any specific 



computation, there is no ambiguity in the physical final results. 

Using the above regularization function, the regularized second order source 
(2,2) 5 z,re g frQm the y = 2 y y = 2 ) coupling is obtained 

= G ^S™% r)6{r - R{t)) + G ( 2 ' 2 )^(t, r)9(R(t) - r) 
+ G ^S s 20 (t,r)6(r - R(t)) + G ™S* (t,r)±6(r - R{t)) 

+ T ^S 5 20 (t, r)S(r - R(t)) + T ^S s 20 (t, r)±-5{r - R(t)) , (61) 

where the superscripts G and T denote the source terms which are derived from 
G^u[h^\ hP^] and T^' h \ respectively. The six source factors, ^'^S^*, etc. are given 
Appendix B| For the factors of the step functions in the above equation, G&tyS^ 



m 



behaves as 0(r 2 ) for large r, and G ( 2 ' 2 ) S\% vanishes as 0((r — 2M) 1 ) at r = 2M. 
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4-3. Regularized second order source from (£ = 0) • (£ = 2) 

The second order source derived from the {£ — 0) • {£ — 2) coupling needs regularization 
of its behavior near the horizon. Note that when we choose another gauge for the first 
order I = mode, some regularization is also needed for the behavior for large r. To 
that end, we use the regularization function, 

X7^\t,r) = 548 s R 1061 r ; 2728 M Hgff (t, r)HST (*, r) 



107 (r-2M)M m d u(1)RW 



2744v^F r iU0U 

+ 274U M(583r > : 756M) g - (V 4 g " r(Vr) ' <62) 
Then, we obtain the second order regularized source from the (£ = 0) ■ (£ = 2) 
coupling 

^Si (t,r) = + (^ (^)-irY 2 f' 2) (^)) 

= G(0 ' 2) 5 °^(t, r)6(r - R(t)) + G(0 ' 2) S^(t, r)6(R(t) - r) 
+ G ^S s 20 (t,r)5(r - R(t)) + G ^S* (t,r)±6(r - R(t)) 

+ T ^S s 20 ( tj r)5(r - R(t)) + T ^S s w (t, r)±5(r - R(t)) 

+ T ^S s 20 (t,r)5(r - R(t)) + T ^S^r)^-8{r - R{t)) . (63) 

Here, we have written the contribution of ha^ a (£ = 2)y o(^)^o*o(^( r )) an d ha^ a (£ = 
0)F 20 (fi)F 2 *(fi(r)) as T{0 ' 2) S 5 20 and T{2 ' 0) ^ , respectively. In the above equation, G (°' 2 )^* 
behaves as 0(r~ 2 ) for large r, and g ^°' 2 ^S 2 q vanishes as 0((r — 2M) 1 ) at r = 2M. 

5. Summary and Discussion 

In this paper, we have completed, for the first time the program of self-consistent binary 
black hole second order perturbations in the small mass ratio limit. Our analysis applies 
to headon collision, but can be extended to arbitrary orbits if worked in the Lorenz 
gauge [51] . The first order perturbations of two black holes starting from rest at a 
finite distance have been solved in the frequency [52] and time [53] domains. Then 
the corrected trajectories (from the background geodesies) via the computation of the 
self- force have been computed in [25l [26] . Here, we obtained the regularized source 
for the second order Zerilli equation in the case of a particle falling radially into a 
Schwarzschild black hole. This is given by fl57|) with flBTj) and (163]) . Using this second 
order source, we are able to compute the second order contribution to gravitational 
radiation by numerical integration of the wave equation (156]) and compare to full 
numerical simulations. The derivation of the second order gravitational wave from the 



second order Zerilli function is discussed in Appendix C 
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A key point is to prove that there is no S 2 term in the second order source. To show 
this, we have used the fact that the first order metric perturbations in the RW gauge are 
C°. In general orbit case (including circular orbits), the first order metric perturbations 
are not C° in the RW gauge, but they are C° in the Lorenz gauge [51]. This later gauge 
choice favors the study of the second order perturbations for generic orbits. 

To be fully consistent in the second perturbative order, we have to include the 
self-force contribution T^ 2 ,' 5 ^ which is derived from the deviation from the background 
geodesic motion. The self-force for a head-on collision has been calculated in [25, 26J, 
and in a circular orbit around a Schwarzschild black hole in [5H [551 ES], but have not 
been obtained in the general case yet. 

Here, we have discussed only the £ = 2 mode of the second perturbative order since 
this gives the leading contribution to gravitational radiation. There remain a question 
about the convergence of the second order metric perturbations, though. Based on the 
works by Rosenthal [571 EEl EH] , we discuss this problem. 

First, we consider a second order wave equation which is given by 

Uh^ = Sf , (64) 

(2) 

where □ is the wave operator and the second order source is derived from the 
first order metric perturbations with a local behavior around the particle location as 
ftM ~ 0(e _1 ) , where the spatial separation e = |x — x | with a particle location x . 
In this discussion, we are not on the particle's world line but rather take a limit to the 

(2) 

particle location. We need second derivatives to compute , then the local behavior 
becomes ~ 0(e~ 4 ) . For the above source, if we solve the wave equation by using a 
usual four dimensional Green's function method, this solution diverge everywhere [57J. 
Therefore, to obtain the finite solution, we need to remove the 0(e -4 ) and 0(e~ 3 ) terms 
from . Here, we note that we can remove the 0(e~ 3 ) terms by using a regular gauge 
transformation [58J. Since we consider the second order gravitational wave at infinity 
which is gauge invariant, the 0(e~ 3 ) terms does not contribute. In practice, we use 
the gauge invariant wave-function in our calculation. Thus, we may discuss only the 
problematic 0(e~ 4 ) terms. 

Here, Rosenthal has already shown the most singular part of the second order metric 
perturbations as 

hf ~ 0(e- 2 ) . (65) 

This is a peculiar solution of \3h^ ~ 0(e -4 ) . Using this solution, we rewrite the second 
order wave solution as 

hf) = - ~ (S^ - 0(e- 4 )) ~ D- 1 (0(e^ 2 )) , (66) 

The right hand side of the last line is finite after the integration by using the retarded 
Green's function as D -1 . The final result is, 

fcOO = hf + h?> , (67) 
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is the physical second-order gravitational perturbations [5B]. Thus, if we construct a 
second order gauge invariant i})^ from the above metric perturbations, this is given by 

= ^ + ^ , (68) 

where ip^ and ip^ are derived from hs and hi 2 \ respectively As a result, the apparent 
divergence derived from the first consideration is only a gauge contribution. 

Next, we consider the expansion in terms of tensor harmonics of the second order 
gauge invariant wave-function (which needs not to be same as the Regge- Wheeler or 
Zerilli function), 

Pfa^S = ^ ■ (69) 
In our situation, we solve the above equation by numerical integration. Formally, we can 

(2) — 1 (2) —1 

write the solution as i/j^ = O im S^ . This \3 em means a numerical integration with an 
appropriate boundary condition. We may also use the retarded Green's function. Since 

(2) 

the solution ip t7 ^ is gauge invariant, this summation over the (£, m) modes coincides 
with ip^ in (168]) . Hence, the summation of ifj^ over modes has a finite value, (except 
for the location of the particle.) In particular, the asymptotic behaviour for large r, 
where we need to compute gravitational radiation, is well defined. 
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Appendix A. The Regge- Wheeler- Zerilli formalism in the time domain 



In this paper, we have discussed the metric perturbations using the Regge- Wheeler- 
Zerilli formalism, with similar notation to that of Zerilli's paper [UJ. There are some 
differences though that we summarize in this appendix. 

For the first and second order metric perturbations, and stress-energy tensors, we 



expand h^l and T) 



1, 2) in tensor harmonics, 



(in 



1 - 



2M 



(m 



1^2l(£+l)h^(t,r)b tm + 



-£(£ +1)(£-1)(£ + 2) 



1/2 



(ill 



+ 



(V2*W(f,r) - 4£+!) G « (t,r)) gem - 



vWTij ,,) 

tl oem\ t i r ) C 0im 
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^iL(*> r ) c &»+ 



[2£(£+l)(£-l)(£ + 2)] 1/2 



2r 2 



Cm 



Im 



(A.l) 



10 



(A.2) 



where a &„, a^ m , • • • are tensor harmonics defined by (3.2-11) in [31]. 

The tensor harmonics can be classified into even and odd parities from the above 
expressions. Even parity modes are defined by the parity (— 1) under the transformation 
(9,4>) — > (tt — 9, 4> + tt), while odd parity modes are by the parity (— Using 
the orthogonality of the above tensor harmonics, we can derive the coefficient of the 
corresponding tensor harmonics expansion. For example, 



A ( oL(t,r) = J T^.a* oem dQ = J 



& m § V ®Tp V Ojoema/3 ^ , 



sin 



(A.3) 

and 8^ a has the component, 



where * denotes the complex conjugate, dVL 
5^ = diag{\, 1, 1/r 2 , l/(r 2 sin 2 9)). 

We summarize the Regge-Wheeler-Zerilli formalism in the time domain. The 
linearized Hilbert-Einstein equation is given in ffTUj) . Next, we specify that the stress- 
energy tensor in the right hand side of the above equation be the one of a point particle 
moving along a geodesic; 



rp(l)lW 



+oo 



\1W 



6®{x- 

DO 

dz» dz u 5{r 



, , . dz>^ dz> T 

- R(t)) 



dt dt r 2 
where we have used the following notations, 



S^{Q-Q(t)) 



# = z»(t) = {T(r), R(r), 9(r), $(r)} , U c 



dT{r) 
dr 



(A.4) 



(A.5) 



for the particle's orbit. This stress-energy tensor is expressed in terms of the tensor 
harmonics as given in table IA1I for the even parity modes. 

Substituting flA.ll) and (1A.2j) into (jTUl) . we obtain the linearized Hilbert-Einstein 
equation for each harmonic mode. Here, we use the RW gauge condition, h 



(i) 

2lm 







for the odd part and h^jj^ = h^j^ = = for the even part. From these ten 



linearized equations, we derive the Regge-Wheeler-Zerilli equations and construct the 
metric perturbations from the Regge-Wheeler-Zerilli functions in the RW gauge. In the 
following, we focus on the even parity modes. 
We introduce the following function, 



2r 



+ 1] 



K 



+2 



in (A r ) 

(r — 2 M) 



(r£ 2 + r£ - 2 r + 6 M) 



H, 



(i) 







2 i'm (A r ) r Kp2 {t, r, 



(A.6) 
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Table Al. The stress-energy tensor for the even parity modes. 



■*£!(*.»■) = 



' ['■§) (r-2M)- 2 5(r-/ 1 '(/iir,;„(()(n) 



= (i- 



(*>»■) = 



0& (*>»■) 



2M\ 



2 <5(r - i?(t))y; m (0(i)) 



^U°—r-H{r - R(t))Y; m (n(t)) 



? m m))/dt 



= [l(£ + l)/^" 1 / 2 ^ f 1 - r -ij( r _ R(t))dY t 

rl /? 

+ l)/2]-VVC/°(r - 2M)- 1 -5(r - R(t))dY e * m (n(t))/dt 
: + !)(*- 1)(* + 2)/2}- 1 / 2 tiU a 5(r - R(t)) 



^V-i2(t)) 



YLm)) 



This function is related to Zerilli's even parity function ipf£ von as 
^(t,r) = rf2 + r£ _ 2 2r + 6M - (r - 2M)^(t.) 

- w,even f , x _ 16 \f2ixi r 2 {r — 2 M) (1) 

The function tp^ n obeys the Zerilli equation, 
d 2 d 2 



dt 2 ^ dr* 2 



4>Z n (t,r) = S^(t,r), (A.8) 



where r* = r + 2Mln(r/2M - 1) and 
r — 9 /If 

™ = r^ P + H-2r + ,Mf + 1)(£ + ^ ' "> V 

+6 M{t + 2) 2 (£ - l) 2 r 2 + 36 M 2 (£ + 2)(£- l)r + 72 M 3 ) , (A.9) 

and the source term is given by 

_ 16 tt (r - 2 M) 2 (rl 2 + rl - 4r + 2 M) (1) 
^ m [t ' r) ~ £(£ + l)( r £ 2 + r£-2r + QM)r ^^r) 

16V2 7r(r-2M)_ y(1) 32 7r(r-2M)V2 ^ 

v/f+lMM ' m (r£ 2 + r£-2r + 6M)v^(7Tl) £m 

32vr(r-2M) 3 x _ ( 16 m j 2 2,3 _ .2,2 

(r^+rf-2r + 6JlO^ + l) 0r U *' rj \ ib7rr ^ r + lr 1 br 1 

+16 r£ 2 M - 6 r 2 £ + 16 r£M + 8 r 2 - 68 rM + 108 M 2 )/[(£ + l)£(r£ 2 + r£ 

. „^2il AD , x 327r(r-2M)r 2 d „ (1) , , 

-2 r + 6 M) } }4L(t, r) + ^ - - -L - - - J— __ - ^ r) 

32y^7r(r- 2 M) 2 (1) 
+ ( r £ 2 + ri _ 2r + 6 M)£{£ + lf £m 1 ' J ' 1 J 



Regular second order perturbations 19 

Here T$> u = have been used to simplify the source term. Using the function ip£^ n , 
the four coefficients for the metric perturbations in the RW gauge are expressed as 

K fl RW (t, r ) = \{ ( £V + 2 r2£3 ~ r2£2 + 6 r ^ M - 2 r 2 £ + 6 rlM 

-12rM + 24 M 2 ) / [(r£ 2 + r£-2r + Q M)r 2 ] } ^ en) (t, r) 

(r-2M) d_ (even) 32^ (1) 

+ r 0r^ m [ ,n £(£+l)(rP + r£-2r + QM)^ oem[ ' J ' 

*£T(t,r) = _ i (g ± ^m- 6 m) ^ (t> r) + r a ^ (t> r) 

1 l(l+l)(rl 2 + rl-2r + 6M) /(even) 
+ 4 ( r -2I)r ^ m lt ' rj ' 

jyt^ff r ) - H {1)RW (t r) + 16 F {1) (t r) 

^(Dw f , r(r-3Af) 9 , __r^ ^ , 

1£m lt ' j " (r-2M)f(l+l) eft £m [ > )+ £(£ + !) dtdr £m l ' j 



Appendix B. Second order source terms 

In this appendix, we show the regularized source for the second order Zerilli equation 
in fl56l) from the contribution of the (£ = 2) ■ (£ = 2) coupling explicitly. We can also 
derive the source from the (£ = 0) ■ (£ = 2) coupling in the same way, but the expression 
is so long that we do not have the space to show it here. The regularized second order 
source (IBT!) from the (£ = 2) ■ (£ = 2) coupling has the following six factors 

3 (8r 2 + 12 Mr + 7M 2 ) (r — 2 M) 



G(2 ' 2) ^*(t,r) = ^V5(r-2M) 



r 2 (2r + 3M) 



d 2 TOU tu ^ 3(2r 2 -M 2 )(r-2M) ( d 2 TOUt . , 



dtdr M y ' 7 r(2r + 3M) 

5 A (18 r 3 — 4 r 2 M — 33 r M 2 — 48 M 3 ) /fl 



X U^' r V " r 2 (2. + 3M) ~ U*^' r; 

x (J^m) "(r-2M) 2 r )) (|-*#(f,r)) - 3(112 r 6 

+480 r 4 M + 692 r 3 M 2 + 762 r 2 M 3 + 441 r M 4 + 144 M 5 ) ^ *ao*(*> r ) 
x r) /(r 3 (2 r + 3 M) 3 ) - feLzMfllzZjO (| r) ) 

52 ,t, ou ^^ . 3(r-2M) 2 M / 9 3 t A 



dtdr 2U v ' 7 r(2r + 3M) V<^<9r 
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d 3 



+(r-2M) 2 —y™(t,r) ^^(t,r) 



dr 3 







dt 



G(2 ' 2) ^(t,r) = -v / 5(r-2M) 
d 



3(r-2M) 2 M 
r(2r + 3M) 

/ d 3 



/[v^(2r + 3M) 
/ <9 3 
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(B.l) 



V dt dr 2 



r-2MY ^—y%(t,r)) \^- 2 r) j - (18 r 3 - 4r 2 M - 33r M 2 

-48M 3 )/(r 2 (2r + 3M)) (jU»(t,r)) r )) 



(r-2M) (3r-7M) / <9 2 



<9r 



#S(t,r) - 3 (112r 5 + 480 r 4 M 



+692 r 3 M 2 + 762 r 2 M 3 + 441 rl 4 + 144 M 5 ) /(r 3 (2 r + 3 M) 3 )^,"(t, r) 



,2 ji j 3 



3 (8r 2 + 12Mr 



r 2 (2 r + 



7M 2 )(r-2M) (JL*^ ^ 
3M) U^r 2o( 7 



3(2r 2 -M 2 ) (r-2M) / 9 



9 3 



r(2r + 3M) 



<9 



<9r 2 



*S5(*,r) 



9* 



/[v^(2r + 3M) r 2 ], 



(B-2) 



G(2 ' 2) ^ (t,r) = -^5^/^(0,0) 



d 3 

x(^3^o(^; 



cLR(t) 



'-R(t) + 2 M) 2 (2M + 3£ 2 i?(t) - #(£)) 



p , , - 2 ^-^(-38 R(t) M 3 — 48 M 4 — 40 i?(t) 2 £ 2 M 2 



+282 M R(t) E + 288 M E + 171 i?(t) 2 M + 66 M E R(t) 3 + 12 R(t) 



d 

-94 R(tfM- 64 R(tfE 2 ) x ( — ^(t,r 



r = i2(t) 



(i?(t) (2i2(i) + 3Mf 



+6 ^^(-48 £ 2 i?(t) 5 + 8 i?(t) 5 + 288 M 5 £ 2 + 546 i?(t) M 4 £ 2 + 52 i?(t) 4 M E 2 
dt 

-60 i?(t) 4 M - 41 i?(t) 2 M 3 + 288 i?(t) 2 M 3 £ 2 + 104 ^(t) 3 M 2 + 36 M 5 
-8 i?(t) 3 M 2 E 2 )¥ 2 n (t, R(t)) j (R(t) 2 (2 R(t) + 3 M) 3 ) - (-R(t) + 2 M) 

x (18 M 3 + 5 i2(t) M 2 + 6 M 2 R(t) E 2 + UM R(t) 2 E 2 - 43 R{t) 2 M + 18 i?(t) 3 

^ d 2 



-12 R(t) 3 E 2 



(<9r dt 



r = R(t) 



(R(t) (2 R(t) + 3 M )) + (216 M 



+162 i?(t) M 4 - 72 i?(t) M 4 E 2 - 150 i?(t) 2 M 3 E 2 + 165 i?(t) 2 M 3 — 4 i?(f) 5 
-386 ^(t) 3 M 2 + 36 i?(t) 3 M 2 E 2 - 108 i?(t) 4 M E 2 + 126 ,R(t) 4 M 



+56E 2 R(t) 5 )(^^(t,r^ 



r = R(t) 



(R(t) 2 (2 R(t) + 3 M) ; 



Q3 



<9r 2 <9t 



r = R(t) 



— (-R(t) + 2 M) 2 (10 M + 6 £ 2 - 5 R(t)) 

+(-R(t) + 2 M) ^^(48 M 3 + 26 R(t) M 2 + 66 M 2 R(t) E 2 + 34 Mi?(t) 2 £ 2 
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13 R(t) 2 M + 12 R(tf E 2 - 6 R{tf) [ ^- #™(f , r 



dr 2 



R(t) 



(R(t) (2R(t) + 3M)) - 



dR(t) 
dt 



-R{t) + 2 Mf (2 M + 3 E 2 R{t) - R(t)) 



Q3 



x 



dr 3 
d 3 
dr 2 dt 



r = R{t) ~ g (-*(*) + 2 M ) 2 (10 M + 6 £ 2 i2(t) - 5 i2(f )) 

r = ^ - 47r ^r 71 y; ° (0 ' 0)(810 m6 + 1113 m5 



+972 M 5 £ 2 + 1710 M 4 i?(t) 2 £ 2 + 1161 M 4 i?(t) 2 + 2124 M 3 R(t) 3 E 2 
-2352 M 3 i?(t) 3 + 328 M 2 R{tf - 1296 M 2 R(tf E 2 - 96 M R(tf E 2 
+120 MR{tf + 32 i?(t) 6 - 288 E 2 R{tf)E 



(R(t) (2 R(t) + 3 M) 4 + 2 M)) 



[£7i2(t) 2 (2i2(*) + 3M)], 



(B.3) 



G ( 2 %^( { ,r) = ^^^(0,0) 



-R(t) +2M)(-10M 2 + 6£ 2 M 2 + 17 R(t)M 



-4 M £ 2 ) - 6 i?(t) 2 + 4 £ 2 it>(t) 2 ) ^ (|- *™(t, , ))r = R{f) 
! {2 R(t) + 3 M) + 6 (-#(*) + 2 M) (3 E 2 M 3 + 2 M 3 + 6 M 2 R(t) E 2 - 3 i2(t) M 2 



+4 M R{tf E 2 + 5 R{tf M + 4 i?(t) 3 E 2 — 2 R(t) 



x 



(i2(f) (2 + 3 M) 2 ) + + 2 M) 2 (2 £ 2 - 3 R(t) + 6 M) 

/ <9 2 



V dr dt 

d 2 



r = R(t) + 2 (_jR(t) + 2 M)2 R(t) " ^ + 2 M) 



i2(f) 



+ (— + 2 M) (3 M 2 + 6 i2(t) M — 2 R(t) 2 ) 



x (2 E 1 i?(t) - 3 R(t) + 6 M) I — tt25(f , r 



5 



9* 



R{t) /(R(t) (2R(t) + 3M)) 



[E R{tf (2 i2(t) + 3 M)] , 



(B.4) 



^fi&ftr) = -VSv^/^ W) -2M)y 2 *(0,0) 



-2(24M 3 + 10 M 2 



+6 M 2 R(t) E 2 + 9M R{t) 2 E 2 — 7 R{t) 2 M + R{tf E 2 — 2 R{t) 3 ' dR ^ 



dt 



d 



dR(t) 



= R{t) I Wt) (2 m + 3 M )) + 6 dt^ 2A M 
+6 M 3 i2(t) E 2 + 36 i2(t) M 3 + 21 R(tf M 2 + 13 R{t) 2 E 2 M 2 + 10 R{tf M 
+8 M E 2 R(tf + 6 i?(t) 4 E 2 - 4 R(t) 4 )V™(t, R(t)) I [R(tf (2 + 3 M) 2 ) 



+ (£ 2 - R(t) + 2M) (-R(t) + 2M) 



d 2 



dr dt 



R(t) 



Regular second order perturbations 



22 



+ (E 2 R(t) - R(t) + 2 M) 



dR(t) 
dt 



6 (E 2 R(t) - R(t) + 2 M) (R(t) 2 + R(t) M + M 2 ) ( d 



R(t) (2R(t) + 3M) 
/{ER{tf{2R{t) + 3M)}, 

T ^S 5 2 ' Q (t, r) = ^^ (R(t) - 2 Mf Y* (0, 0) 



r = i?(t) 



(B.5) 



dR(t) 
dt 



'-R(t)+2M) 



_ rf.R(t) 6 (#(t) 2 + M + M 2 
r = R(t) dt R(t) (2R(t) + 3M) 



**S(t, i2(t)) 



(E R(t) 3 (2 R(t) + 3M)) . 



(B.6) 



Here, we have used ( 127)) for \l/f and its derivatives at the particle trajectory. And for the 
homogeneous solutions, ^o* an d ^20 have been used to write the above source terms. 

Appendix C. Second order gauge transformation 

In this appendix, we deal with first and second order gauge transformations. In order 
two obtain the second order waveform, it is necessary to derive the second order metric 
perturbations in an asymptotic flat (AF) gauge. The Regge-Wheeler-Zerilli formalism 
that we have employed in the RW gauge is not asymptotically flat. Therefore, we will 
focus on the gauge transformation from the RW gauge to an AF gauge. We also need 
to discuss the first order gauge transformation to an AF gauge simultaneously. 
Here, we consider the following gauge transformation [351 H2J . 



x 



RW 



1" 



AF 



X 



,H -I" ] + , (C.l) 



where comma "," in the index indicates the partial derivative with respect to the 
background coordinates, and fp 1 ^ and £( 2 )^ are generators of the first and second order 
gauge transformations, respectively. Then, the metric perturbations changes as 
1.(1) 1.(1) _ 1.(1) __ r 



(C.2) 
(C.3) 



1.(2) 1.(2) _ t(2) 1 r \^r 2 /,(!) 

Next, we discuss the (£ — 2) ■ (t — 2) and (t — 0) • (£ — 2) parts separately. 

Appendix C.l. First order i = 2 mode and second order (£ = 2) • (I = 2) part 

In this paper, we have used only the even parity mode, therefore a generator of the 
gauge transformation for I = 2, m = modes can be written as 

(t, r)Y 20 (9, 0), Vf (t, r)Y 20 (9, 0), 



S£=2 



Vf (t,r)«9 e F 20 (M), Vr(t,r) 



(0/ 



sin 2 6 



(C.4) 
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where i = 1 and 2 denote the first and second perturbative order, respectively. There 
are three degrees of gauge freedom for each perturbative order. 

The gauge transformation of the metric perturbations is explicitly given as follows: 
For the first order metric perturbations, we find 

H®T(t,r) = ff«r («, r) + 2 JV„« (t, r) + 2 r(r f 2M) Vl"» («, r) , 



rj 

r <9r v ' 7 r — 2 M dt 



H^ F (t,r) = FST(*,r) -2^) (t,r) + 2^^^^ (t 

= (t,r) - r 2 |v 2 « (f,r) , 

h[T AF (t,r) = -j^VP (t,r) - r^V^ (t,r) , 

G? )AF (t,r) = -2y 2 «(t,r) . (C.5) 

For the second order metric perturbations, we can calculate the gauge transformation 
straightforwardly, but we obtain very long expressions. For example, they can be written 
formally as 

KT F (t, r) = KT W (t, r) - V> (*, r) + 5K^(t, r) , (C.6) 

r 

h^ )AF (t,r) = -^jj^V^ (t,r) - r ^V 2 ^ (t,r) + 5h^\t,r) . (C.7) 

GT F (t, r) = -K 2 (1) (t, r) + tfGg(t, r) , (C.8) 
where 8K^ , Sh^l^ and 5G^q are defined by the tensor harmonics expansion of the 



last two terms in the right hand side of ( 1C.3I) . This includes only quadratic terms of the 
first order wave-function. 

First, we consider the asymptotic behavior on the i = 2 mode of the first order 
metric perturbations in the RW gauge. The asymptotic expansion of the wave-function 
ip2o en is given by 

+0(r- 3 ) , (C.9) 

where we have introduced T r = t — r*(r) for simplicity. In the following calculation, we 
need only the leading order contribution with respect to the above large r expansion. 
Then, the coefficients of the metric perturbations are given, from ( lA.llj) . as follows 

H$r(t,r) = H^T (i,r) = ± (J^w) r + 0(r°) , 



V^^r) = -- (-^F(T r ))r + 0(r°), 
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KT W (t,r) = -\ Jjn>F(T r ) + 0(r- 1 ) . (CIO) 

On the other hand, the metric perturbations in an AF gauge should behave as 

4T(t, r) = H$ AF (t, r) = h^ AF (t, r) = . H^t, r) = 0(r~*) , 

h[T AF (t,r) = 0(r^),K$ AF (t,r) = O(r^) , G$ AF (t,r) = O^ 1 ) . (C.ll) 

This asymptotic behavior will also be considered for the second order calculation. We 
find the following gauge transformation to go to the AF gauge. 

1 / d 3 

6 \dT r * 

Vl 1 \t,r) = - 1 - (J^F (r r ))r + 0(r°) , 

V} 1 \t,r) = - 1 - (j^F(T r )y-i + 0(r- 2 ). (C.12) 

The above results are calculated iteratively for large r expansion. Since the transverse- 
traceless tensor harmonics for the even parity part is fi m in ( lA.lj) . the coefficient of the 
metric perturbations related to the gravitational wave is G < f^ F . This becomes 

GT F (t,r) = ll^F(T r ) + 0(r->) 

= ^ 2 e nt,r) + 0(r- 2 ). (C.13) 

with the use of flClQl) 

Next, we discuss the second perturbative order. When we treat the second order 
metric perturbations from the (£ = 2) • (£ — 2) coupling, in practice, we calculate xfo 
numerically instead of x| (t,r), where xfo nas been considered in (1551) as 

Xl(t,r) = x 2 z (t,r) - X r 2 o g ' (2 ' 2) (t,r) - &™(t,r) , (C.14) 

where X20 * s the (I — 0) • (£ — 2) contribution, to be discussed in the next subsection. 
Here, to derive the gravitational wave amplitude for the second perturbative order, we 
also need to obtain the coefficient G^ AF in an AF gauge as in the first order case. 
The asymptotic expansion of xfo * s 

Xlofrr) = ± ^F 2 (T r ) + (Xr- 1 ) . (C.15) 

Here, we may consider only the leading order contribution with respect to large r 
expansion in the same manner as the first order calculation. The xio contribution 
to the waveform is derived by the same method as that for the first perturbative order. 
First, we obtain dK^ RW /dt in the RW gauge from (}5Tj) as 

d „(2)rw /, x 1 d A „ , m s , y/5 ( d 4 „ , m N \ d 3 



t-K^' (t, r) = —^F 2 (TV) + — jF ( TV) — — ( TV) 

dt 20 v ' ; 3 dT r 4 v ; 18v^F \dT r A v 7 dT 3 v ; 

+0(r" 1 ), (C.16) 
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where the second term in the right hand side of the above equation arises from the 
regularization function Xsjcf'^ anc ^ the anc ^ £>2(? "terms i n (1511) ■ Integrating the 
above equation for W 



K. 



(2)RW 
20 



(t,r) 



1 d 3 



3 dT r 



■F 2 (T r ) 



V5 



d 3 



r- 1 -^f(t; 



+ 0(r 



-1> 



(2) 

The second order gauge transformation of K^q is given by 



K. 



(2)AF 
20 



(t, r ) = Kg* W (t, r) - V } (t, r ) + (t , r) 



where, (Sii^o is defined by the tensor harmonics expansion of (l/2)£ 2 w g fiu — JC^h 1 ^-^ 
in (IC.3I) and derived as 

2" 



(C.17) 



(C.18) 



^ 2) (t,r) 



d 1 



dT r 



rF(T r 



d 2 



dT r 



>F ( T r ] 



d 3 
~df} 



F ( T r ) 



(C.19) 



2520F 
+0(r- 1 ). 

From the above results and the AF gauge condition for K^j in (IC.11I) . in order to remove 
the 0(r°) terms, the second order gauge transformation Vi 2 ^ is 



1 

3 



dT r 



:Fo T r 



V5 



252^ 



d 3 
dTr 



F (Tr 



+ 



d: 1 



dT r 



- F ( Tr 



d 2 

dT r 2 



r + 0(r c 



(C.20) 



Here, we note that it is sufficient to consider the leading order with respect to large r 
to derive the second order waveform. 

Next, we consider to derive V^ 2 -* from the condition of ^120^- The second order 
gauge transformation is given by 



h (e)(2) AF 
"'120 



(t,r) 



(C21) 



dhxl^ is defined by the tensor harmonics expansion of (l/2)£ 2 w g fll/ — C^)h^ Wilv in 
Eq. (1(130 . By considering the asymptotic expansion, we obtain 



Sh^\t,r) 



%/5 



d A 



dT r 



{F {T r ) 



d 2 



dT r 



jF(T r ) + 



d 3 
dTr 



:F(T r 



5O40F 

+0(r°). (C.22) 

Then, V^ 2 -* is calculated from the above value and the result for in (IC.20I) with the 
AF gauge condition in Eq. (jC.llh as 



|_K 2 (2) ( t ,r) = --^F 2 (T r ) r- 1 
dr v ' ; 3 dT r 3 V ; 

d 



V5 



d 3 



-*F(Tr 

630F \dT r 3 



r 1 + 0(r 



-2\ 



= -^2 (2) (t,r) + 0(r- 2 ). 
In the last line, we have used the definition of the retarded time, T r — t — r*(r). 



(C.23) 
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From the above results, we can consider the metric perturbations related to the 
gravitational wave amplitude, i.e., F ■ The gauge transformation is given by 

GT F (t,r) = -V 2 ^ (t,r) + 6G$(t,r) . (C.24) 



Here, SG^j is defined by the tensor harmonics expansion of (l/2)£ 2 (1) g tlL/ — C^)h!^ w ^ v 
in Eq. (IC.3I) . and found to be 

y/b ( d 2 



rF(T r 



6G®(t,r) = —— j F ( T r ) , , . 

20 v ' ; 126^ \dT r 2 K 1 ) \dT 3 

Inserting fl023|) and (10251) into (l024j) . we obtain 



d 3 



o( r - 



(C.25) 







1 d 3 



dt G " o(t ^ 3dT r 

V$ \ ( d 3 



3^ mr- 1 



126v^F 



dT r 



■F (T r 



+ 



d: 1 



dT r 



dT r 



rF(T r ) 



- 1 + 0(r~ 2 ). (C.26) 



The gravitational waveform is now (by use of ip2o en an d X20 
^G%{t,r) = ~xl{t,r) 



V5 1 



+ 0(r" 2 ) . 



(C.27) 



14a/7t r 

If we consider higher order corrections with respect to the 1/r expansion, we can show 
that all metric components satisfy the asymptotic flat gauge condition in fIC.llj) . 



Appendix C.2. Second order (£ = 0) • (£ = 2) part 



We have already discussed the first order perturbations for the £ = mode in section [3^41 
In this paper, we use the first order £ = metric perturbation given by (jUJ). This 
satisfies the AF gauge condition in (IC.llj) . Therefore, it is not necessary to consider 
the first order gauge transformation of the £ = mode in flC.2j) and (IC.3I) . and we will 
focus on the second perturbative order related to the (£ = 0) • (£ — 2) coupling. 

We discuss the gravitational wave amplitude for the second perturbative order 
which arises from the (£ = 0) ■ (£ = 2) coupling by using the same method as in the case 
of the (£ — 2) • {£ — 2) part. Note that we have already discussed the contribution from 
xfo an d X20 ,< " 2 ' 2 ' > i n the second order wave function of ( |55l) . In the following, we consider 
the 

X2cf an d the contribution of the last term in the right hand side of (1C.3I). i.e., 
— ^mhftw where ^ is the generator of the gauge transformations for £ = 2, and we 



(i) 

/eg,(0,2) . 



use (H4l) as the £ = mode of h 

The regularization function x-zcf '^'^ ^ s gi yen i n fl62|) . This function behaves 0(r" 
for large r, therefore we expect that X20 '^°' 2 ^ does not contribute to the second order 
gravitational waveform at infinity. 

In the same way as for the {£ — 2) • (i — 2) part, we consider the gauge 

Here, we obtain 



transformation of in ( 1C.18 



0{r 



-3\ 



(C.28) 
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where K^j RW and SK^ arise from Xsacf'^ anc ^ the tensor harmonics expansion of 
— C'^mh^r in ( 1C.3I) . respectively. This is already asymptotically flat, therefore we do 
not need any further gauge transformation, for the (£ = 0) ■ (£ = 2) part. 

^( 2 )(t,r) =0(r- 2 ). (C.29) 

In (E2U), we have 

6h^ e \t,r) = 0( r - 3 ), (C.30) 

where Shfl^ is defined by the tensor harmonics expansion of — C^h^y in ( 1C.3I) . 
From this equation, we also conclude 

l.V 2 V(t, r )=0(r- 4 ). (C.31) 

Hence, there is no contribution from the [I = 0)-(£ = 2) coupling to the second order 
gravitational wave, except for x| . This means that if we obtain x| in the numerical 
calculation, we can obtain the gravitational wave amplitude for the second perturbative 
order by using (1C.27I) . 
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